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Qh' Abstract 

T^, - The three-dimensional XY model is studied at finite chemical potential 

using complex Langevin dynamics. The validity of the approach is probed 
\^J I at small chemical potential using imaginary chemical potential and conti- 

QQ ■ nuity arguments, and at larger chemical potential by comparison with the 

^O . world line method. While complex Langevin works for larger /3, we find 

that it fails for smaller f3, in the region of the phase diagram corresponding 
to the disordered phase. Diagnostic tests are developed to identify symp- 
(<— ^ ■ toms correlated with incorrect convergence. We argue that the erroneous 

^^ i behaviour at smaller /3 is not due to the sign problem, but rather resembles 

dynamics observed in complex Langevin simulations of simple models with 
complex noise. 
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1 Introduction 

Field theories with a complex action are difficult to treat nonperturbatively, be- 
cause the weight e~^ = \e~^\e^'^ in the partition function is not real. Standard 
numerical approaches based on a probability interpretation and importance sam- 
pling will then typically break down, which is commonly referred to as the sign 
problem. This is a pressing problem for QCD at nonzero baryon chemical poten- 
tial, where a nonperturbative determination of the phase diagram in the plane 
of temperature and chemical potential is still lacking [T]. Several methods have 
been developed to explore at least part of the phase diagram p]-[T2]. but in gen- 
eral these can only be applied in a limited region. Recent years have also seen 
an intense study of the sign problem in QCD and related theories, which has led 
to new formulations [13l[14] and considerable insight into how the complexity of 
the weight interplays with physical observables [T5] - [2^ . Finally, in some theo- 
ries the sign problem can be eliminated completely, using a reformulation which 
yields a manifestly real and positive weight [251428] . This demonstrates that the 
sign problem is not a problem of principle for a theory, but instead tied to the 
formulation and/or algorithm. For QCD an exact reformulation without a sign 
problem has unfortunately not (yet) been found. 

Complex Langevin dynamics [2^150] offers the possibility of a general solution 
to this problem. In this formulation the fields, denoted here collectively as </>, are 
supplemented with an additional fictional Langevin time 'd and the system evolves 
according to the stochastic equation. 



In the case that the action is complex, the fields are complexified as 

</)^</.^ + z0\ (1.2) 

and the Langevin equations read (using general complex noise) 

^'^^ Kf+./l^V^, Kf = -Re^_^^, (1.3a) 



<^->-0f'-+i0^ 






The strength of the noise in the real and imaginary components of the Langevin 
equation is constrained via N-^ — Ni = 1, and the noise furthermore satisfies 

iv^m = {vim = ivfi^Hm = 0, (1.4a) 



i.e., it is Gaussian. Since the complex action is only used to generate the drift 
terms but not for importance sampling, complex Langevin dynamics can poten- 
tially avoid the sign problems 

In the limit of infinite Langevin time, noise averages of observables should 
equal the standard quantum expectation values. For a real action/Langevin 
dynamics, formal proofs that observables converge to the correct value can be 
formulated, using properties of the associated Fokker-Planck equation [37]. If 
the action is complex and the Langevin dynamics extends into the expanded 
complexified space, these proofs no longer hold. Nevertheless, a formal derivation 
of the validity of the approach can still be given, employing holomorphicity and 
the Cauchy-Riemann equations. We sketch here the basic notion, suppressing all 
indices for notational simplicity, and refer to Ref. [15] for details. 

Associated with the Langevin process (11. 3p is a (real and positive) probability 
density P[(f)^,(j)^;'d], which evolves according the Fokker-Planck equation 
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with the Fokker-Planck operator 
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Stationary solutions of this Fokker-Planck equation are only known in very special 
cases |33p43 [ 145 ( 148]. Expectation values obtained by solving the stochastic process 
should then equal 
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However, we may also consider expectation values with respect to a complex 
weight p[0; ■(?], 

where, using Eq. (II. ip . p evolves according to a complex Fokker-Planck equation 
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This equation has the desired stationary solution p[0] ~ exp(— S"). 

Under some assumptions and relying on holomorphicity and partial integra- 
tion [l5], one can show that these expectation values are equal, and 
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^Early studies of complex Langevin dynamics can be fomid in, e.g., Refs 
contains a further guide to the hterature. More recent work includes Refs. [38H4^ 



(1.10) 
Ref. [37] 



If it can subsequently be shown that 



hm (O)p(tf) = (O)p(oo), p(0; oo) ~ exp(-S'), (1.11) 

W— 7>00 



the apphcabihty of complex Langevin dynamics is demonstrated. In Ref. |15] this 
proposal was studied in some detail in the case of simple models. Remarkably 
it was found that for complex noise (A^i > 0), the Langevin dynamics does not 
converge to the correct answer. On the other hand, for real noise (A^i = 0) correct 
convergence was observed. 

In this paper, we continue our investigation into the applicability of complex 
Langevin dynamics at finite chemical potential pT| - H5] . We consider the three- 
dimensional XY model for a number of reasons. We found earlier that this theory 
is very sensitive to instabilities and runaways and therefore requires the use of 
an adaptive stepsize [H]. This is similar to the case of QCD in the heavy dense 
limit [mm]. As QCD, this theory has a Roberge- Weiss periodicity at imaginary 
chemical potential |^H9]. Furthermore, it is closely related to the relativistic 
Bose gas at finite chemical potential, for which complex Langevin dynamics was 
shown to work very well (at weak coupling in four dimensions) |l2lll3]. Finally, 
this theory can be rewritten using a world line formulation without a sign problem 
[271128] . which can be solved efficiently using the worm algorithm [28l[50]. This 
allows for a direct comparison for all parameter values. 

The paper is organized as follows. In Sec. |2l we remind the reader of some 
details of the XY model at real and imaginary chemical potential, the adaptive 
stepsize algorithm we use and the related phase- quenched XY model. The world 
line formulation and some properties of the strong-coupling expansion are briefly 
mentioned in Sec. [31 We then test the validity of complex Langevin dynamics in 
Sec. HI and develop diagnostic tests in Sec. [51 In the Conclusion we summarize 
our findings and discuss possible directions for the future. 

2 XY model 

The action of the XY model at finite chemical potential is 

2 

5' = -/3^^cos(03; -0^+p-2/i5^,o), (2.1) 

where < (px < 27i. The theory is defined on a lattice of size Q = NrN^, and 
we use periodic boundary conditions. The chemical potential // is coupled to 
the Noether charge associated with the global symmetry (p^ -^ (px + C( and is 
introduced in the standard way [51]. The action satisfies S*{fi) = S{—fi*). At 
vanishing chemical potential the theory is known to undergo a phase transition 
at /3c = 0.45421 [25l[5^ between a disordered phase when /3 < /3c and an ordered 
phase when /3 > /3c. 



The drift terms appearing in the complex Langevin equations are given by 
^x = -P Yl[^^'^^^^ ~ ^^+^^ cosh(0^ - 0^+^ - /i5^,o) 

(2.2a) 





+ sin(0^ - </.^_,) cosh(0^ - <Pl_, + /i5,,o) 
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-/3 > ; cos(0^ - 0^+,) sinh(0^ - 0i+, - /i5,,o) 




+ cos(0^ - 0^,) sinh(0^ - 0i_^ + /i5,,o) 



(2.2b) 

The equations are integrated numerically by discretizing Langevin time a.s "d = 
nen with e„ the adaptive stepsize. Explicitly, 

0^(n + 1) = 0^(n) + enKfin) + ^Vxin), (2.3a) 

<l>lin + l)=<l)lin) + e^Klin), (2.3b) 

where we specialized to real noise, with {rix{n)r]x'{n')) = 25xx'5nn'- In the case 
that /i = 0^ = 0, the drift terms are bounded and \K'^\ < 6(3. When 0^ ^ 0, the 
drift terms are unbounded, which can result in instabilities and runaways. In this 
particular theory, much care is required to numerically integrate the dynamics in 
a stable manner and we found that an adaptive stepsize is mandatory [H]. At 
each timestep, the stepsize is determined according to 

'n = ^^[^,e^j^j, (2.4) 

where 

ir^"" = max \Kf{n) + iKl{n)\ . (2.5) 

X ' ' 

Here e is the desired target stepsize and {K™^^) is either precomputed or com- 
puted during the thermalisation phase. All observables are analyzed over equal 
periods of Langevin time to ensure correct statistical significance. 

The observable we focus on primarily in this study is the action density {S) /Vt. 
After complexification the action is written as 5* = S^ + zS'\ with 



s""- 




(2.6a) 


s'- 


= (3 Y^ sin(0^ - 0^+,) sinh(0^ - 0^+, - /i5.,o). 


(2.6b) 



After noise averaging, the expectation value of the imaginary part is consistent 
with zero while the expectation value of the real part is even in /z, as is expected 
from symmetry considerations. 

By choosing an imaginary chemical potential // = i/ii the action (12. ip becomes 
purely real. This has both the advantage of enabling standard Monte Carlo 



algorithms to be applied (we choose to employ real Langevin dynamics) and that 
the behaviour at /i^ > can be assessed by continuation from /i^ < 0. The action 
and drift term with imaginary chemical potential are 

5'imag = - /5 X^ COS(0^ - (p^+o + fJ,l5^fi), (2.7) 

Kx = - (3^^ [sin(0^ - (f)^+0 + fii5„fi) + sin(0^ - (f)^_o - l^iS^fi)] . (2.8) 
1/ 

This theory is periodic under /ij — )■ /ij + 271 /N^-, which yields a Roberge- Weiss 
transition at /ij = tt/Nt-, similar as in QCD [15] • This periodicity can be made 
explicit by shifting the chemical potential to the final time slice, via the field 
redefinition 0x,r — >■ 4>xt — 'Px,t — /^iT- The action is then (for arbitrary complex 
chemical potential) 

5'fts = -/5 5Z ^°^('^^' ~ '^^•+'> ~ iNr^i'^T.Nr^ufl)- (2.9) 

x,v 

We have also carried out simulations with this action and confirmed the results 
obtained with the original formulation. The sole exception was the largest /3 
value (/9 = 0.7), where the original action missed the Roberge- Weiss transition, 
while the final-time-slice formulation located it without problems. 

The severity of the sign problem is conventionally (see e.g. Ref. [1]) estimated 
by the expectation value of the phase factor e*"^ = e~^ /\e~^\ in the phase quenched 
theory, i.e. in the theory where only the real part of the action (12. ip is included in 
the Boltzmann weight. In this case, the phase quenched theory is the anisotropic 
XY model, with the action 

Spq = - ^ Py cos{(j)^ - (j)x+o), (2.10) 

x,u 

where (3q = /3cosh/i, and (3i^2 = P- 

3 World line formulation 

The advantage of the XY model is that it can be formulated without a sign 
problem by an exact rewriting of the partition function in terms of world lines 
[271 [28] o Moreover, this dual formulation can be simulated efficiently with a 
worm algorithm [2511^]] . which allows us to compare the results obtained with 
complex Langevin dynamics with those from the world line approach. We briefly 



^The world line formulation has of course a long history in lattice gauge theory, see e.g. 
Ref. [53]. Recent work includes Refs. [54ll56) . For a review, see Ref. [27] , 



repeat some essential elements of the world line formulation and refer to Ref. [28] 
for more details. The partition function can be rewritten using the identity 



g/3 cos 



E W)e""^, (3.1) 



where /fc(/3) are the modified Bessel functions of the first kind. Using this re- 
placement and integrating over the fields, the partition function is written as 

Z = I Dcl^e-' = Y.X{^^^M^'^'"''-''''^ ( E [^-''^ - ^^~^A ■ (3.2) 

•^ \k] x,u \ u J 



The sum over [k] indicates a sum over all possible world line configurations. Since 

dp 



(S) = —(^^^§^-i the action can be computed from 



(^) = -/^(E 



4.,.-i(/3) K 



KM P 



(3.3) 

wl 



where the brackets denote the average over world line configurations. To com- 
pute this average, we have implemented the worm algorithm, following Ref. [2B] . 
We note here, amusingly, that the world line formulation has a sign problem at 
imaginary chemical potential. 

Inspired by Ref. [57], we have also studied a (low-order) strong-coupling ex- 
pansion of this model, using 



k / 2 4 



At strong coupling the chemical potential cancels in most world lines, except 
when the world line wraps around the temporal direction. At leading order in 
the strong-coupling expansion, it then appears in the combination (|/3e'^)^^. In 
the thermodynamic limit it therefore contributes only when |/3e^ > 1. Hence a 
simple strong-coupling estimate for the critical coupling at nonzero /i is given by 

/3,(/i) = 2e"^ (3.5) 

The /i-independence at small /3 and // is known as the Silver Blaze feature in 

QCD [sg. 

The partition function is expressed in terms of the free energy density / as 
Z = exp(— f2/). A strong-coupling expansion to order (3'^ on a lattice with A',- > 4 
yields 

f = -ll3'-f^l3' + Oin (3.6) 



and hence 

(S)/fi = -^/32-^/3^ + 0(/3«). (3.7) 

I ib 

In the phase quenched theory we find 

/pq = -^/32 (2 + cosh^ ii) - -^/3^ (14 + 8 cosh^ /i - cosh^ /i) + 0{p>^). (3.8) 

We can now estimate the severeness of the sign problem at strong coupling. The 
average phase factor takes the standard form, 

(e^^)pq = A = exp [-fiA/] , A/ = / - /p^, (3.9) 

where in this case 

A/ = ^/32 (cosh^ /i - 1) + 7^/3^ (cosh^ /i - 1) (7 - cosh^ p) + C(/36). (3.10) 

On a finite lattice and for small chemical potential we find therefore the sign prob- 
lem to be mild in the strong-coupling limit, since the volume factor is balanced 
by /3VV4 < 1. 

4 Comparison 

We start to assess the applicability of complex Langevin dynamics for this model 
at small chemical potential. In this case we can use continuity arguments to 
compare observables at real and imaginary chemical potential. In Fig. [T]the real 
part of the action density is shown as a function of /i^, for several values of /3: from 
the ordered phase at large /3 to the disordered phase at low /3. We observe that 
at the highest values of /3 this observable is continuous across /i^ = 0, which is a 
good indication that complex Langevin dynamics works well in this region. The 
cusp at /ii = tt/A^t- (corresponding to yU^ = —0.154) refiects the Roberge- Weiss 
transition. At lower /3, however, we observe that the action density is no longer 
continuous: this is interpreted as a breakdown of complex Langevin dynamics. 
In order to verify this. Fig. [1] also contains the expectation values of the action 
density found using the worm algorithm in the world line formalism for real [i. 
As expected, in this case the action density is continuous across //^ = for all 
values of /3, confirming the interpretation given above. We have verified that the 
jump in the action density at lower /3 is independent of the lattice volume. We 
have also verified that the discrepancy at /x^ = between real Langevin dynamics 
and the world line result (at e.g. /3 = 0.4) is due to the finite Langevin stepsize. 
For small /3, the numerical results found with the worm algorithm are con- 
sistent with those derived analytically in the strong-coupling limit above. The 
expectation value of the action density is /i independent and hence the Roberge- 
Weiss periodicity is smoothly realized. Using Eq. (13. 7p . we also find quantitative 
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Figure 1: Real part of action density {S)/Q as a function of //^ on a lattice of 
size 8'^, using complex Langevin dynamics and the world line formulation at real 
fi (/i^ > 0) and real Langevin dynamics at imaginary /i (//^ < 0). The vertical 
lines on the left indicate the Roberge-Weiss transitions at fii = vr/S. 
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Figure 2: Colour plot indicating the relative difference AS* between the expecta- 
tion value of the action density obtained with complex Langevin dynamics and 
in the world line formulation, see Eq. f l4.2p . Also shown is the phase boundary 
/3c (/i) between the ordered (large (3) and disordered (small /3) phase 



agreement: in the strong-coupling expansion {S)/Q = —0.0621 + (9(10 '') for 
13 = 0.2 and -0.145 + O{10-^) for /3 = 0.3. 

As discussed above, for the parameter values and lattice sizes used here the 
sign problem is not severe: taking /i^ = 0.1 and /3 = 0.2, we find that 



2, ,2 



OA/^fi 



/3V^ 



0.51, 



o^'fiX 



/pq 



0.60. 



(4.1) 



We take this as a first indication that the observed breakdown is not due to 
the presence of the sign problem, especially since complex Langevin dynamics 
has been demonstrated to work well in other models where the sign problem is 
severe [121113]. 

To probe the reliability of complex Langevin dynamics for larger values of 
/i, we have computed the action density for a large number of parameter values 
in the (3 — fi plane. Our findings are summarized in Fig. [21 where we show the 
relative difference between the action densities obtained with complex Langevin 
(cl) and in the world line formulation (wl), according to 



AS 



{S)wl — {S)cl 



(S) 



wl 



(4.2) 
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Also shown in this figure is the phase transition fine Pdf^), taken from Ref. |28j . 
We observe a clear correlation between the breakdown of Langevin dynamics 
and the phase boundary: complex Langevin dynamics works fine well inside 
the ordered phase, but breaks down in the boundary region and the disordered 
phase. The largest deviation around /i = 2 is due to the Silver Blaze effect: the 
difference between the action density found with complex Langevin dynamics and 
the correct /i-independent action density is maximal just before crossing over to 
the other phase, where the agreement improves quickly. 

5 Diagnostics 

In this section we attempt to characterize the results presented above in terms 
of properties of complex Langevin dynamics and the distribution P[(f)^, 4>^] in the 
complexified field space, see Eq. f lL7p . We suppress Langevin time dependence, 
since we always consider the quasi- stationary regime, i.e. the initial part of the 
evolution is discarded (we considered Langevin times up to -i? ~ 2 x 10^). Our aim 
is to argue that the discrepancy at small (3 is introduced by complex Langevin 
dynamics rather than by the presence of a chemical potential and hence not due 
to the sign problem. 

A first test of the validity of complex Langevin dynamics is to compare sim- 
ulations at /i = using a cold start, i.e. with (j)^ = initially, and a hot start 
in which 0^ is taken from a Gaussian distribution|j When // = 0, a cold start 
corresponds to real Langevin dynamics. In the case of a hot start, however, the 
fields lie immediately in the complexified space and so the dynamics is complexi- 
fied. Comparison of results obtained with these two initial ensembles gives insight 
into the inner workings of complex Langevin dynamics. We have computed the 
expectation value of the action density at // = using both a hot and a cold 
start. We found them to agree at large /3, despite the fact that the imaginary 
components of the field are initialised randomly. However, when /3 < 0.5, they 
disagree. Moreover, the result from the cold start agrees with the one obtained in 
the world line formulation. One is therefore led to conclude that when /i = the 
imaginary components (f)^ are driven to zero (more precisely, to a constant value) 
at large /3 but are not constrained at small /3. In other words the drift terms 
are not capable of restoring the reality of the dynamics. It is tempting to relate 
this to being in (or close to) the disordered phase. We note that it cannot be 
understood from the classical fixed point structure, since this is independent of /3. 
We also remark that the dynamics at small /3 resembles Langevin dynamics with 
complex noise (A'l > 0) |15], where the trajectories are kept in the complexified 
field space by the stochastic kicks on 0^ (rather than by the drift terms, as is the 
case here). 



•^The real components (j)^ are taken from a Gaussian distribution always. 
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Figure 3: Width of the distribution P[(p^, 0^] in the imaginary direction for var- 
ious values of /3 as a function of /i^ on a 10^ lattice (left) and, for larger fi, as a 
function of /x on a 8^ lattice (right). 
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In terms of the distribution P[(f) ,(j)], these findings imply that P[ 
e~^6{(j)^) at large /3, but not at small /3. This can be further investigated by 
studying the width of the distribution in the imaginary direction^ 



(V)^ 




(5.1) 



When /i = the width should vanish, while when turning on /i one may expect 
it to increase smoothly. The results are shown in Fig. [31 For the larger /3 values 
this is exactly what is observed: the width increases smoothly from zero. For 
the smaller /3 values, however, we observe that the width is nonzero even when 
/i = (when using a hot start), and remains large for nonzero fi. At larger values 
of /i the width is driven again towards zero and agreement with the world line 
results improves, see Fig. [2j We remark here that it is possible that different 
distributions (with different widths) yield the same result for observables. This is 
what is theoretically expected in the presence of complex noise (A^i > 0) [IS] and 
can be seen analytically in gaussian models with complex noise, where a continu- 
ous family of distributions P[0^, 0^; Ni] all yield the same result for observables, 
independent of A^i, even though the width of these distributions is nonzero and 
increases with Ni [59]. In the case we study here, however, we find that the 
failure of complex Langevin dynamics in the disordered phase is correlated with 
the spread of the distribution P[0^, 0^] in the noncompact direction. We con- 
clude that a relatively narrow distribution, with a smoothly increasing width, is 
required. We note again that this resembles observations made in simulations of 
nongaussian models with complex noise [^|60]. 



The mean value 
of initial conditions. 



0; in the large /3 phase, this requires averaging over a large number 
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To investigate the interplay between (the width of) the distribution and ob- 
servables, we express expectation values as 

(A[0^,0i]) = iy'D0^D0ip[0^,0^]A[0^,0i], (5.2) 

with 

Z= [ D(t)^D(t)^P[(t)^,(j)^]. (5.3) 



In general the operator A is not required to be holomorphic, i.e. a function of 
0^ + i(/)\ since this will allow more insight in properties of the distribution|f| The 
distribution of an operator A can then be defined according to 

(A) = I dAPiA)A= ^ I D<p''D<p'P[<p'',<p']A[<p'',<Pl (5.4) 

PiA) = iy"D0«Z}0Ip[0^0I]5(A- A[0^0I]), (5.5) 

with the normalization 

^dAP{A) = l. (5.6) 



where 



Distributions P{A) can be constructed numerically, by sampling A from config- 
urations generated by complex Langevin dynamics. 

The distribution for the action density is shown in Fig. HJ comparing again a 
hot and cold start at /i = 0. This figure supports the earlier claim that real and 
complex Langevin match at larger /3 but fail at smaller /3. However, the reason 
for failure is somewhat subtle. Naively, one might expect a large "tail" caused 
by excursions in the complexified field space to affect the expectation value but 
this does not appear to happen. Instead we find that the entire distribution is 
shifted and becomes only slightly wider at /3 < 0.5 when the hot start is used. 

Finally, the observed difference at large and small /3 also appears prominently 
in the actual dynamics, i.e. in the drift terms. We have analyzed the maximal 
force K"^^^ appearing in the adaptive stepsize algorithm. In the case of real 
Langevin dynamics, the drift terms are limited by an upper bound of K'^^^ < 6/3. 
In the complexified space there is no upper limit and the drift terms can in 
principle become several orders of magnitude larger |13]. The distribution of K^^^ 
is plotted at /i = with hot and cold starts in Figure [51 In the large (3 phase, 
the distributions appear identical, with K^^^ < 6/3. This is consistent with the 
conclusion reached above. In the low /3 phase the distributions are dramatically 
different: in the complexified dynamics, triggered by the hot start, much larger 
forces appear. The distributions are no longer peaked but very broad with a long 
tail (note the horizontal logarithmic scale). At /3 = 0.5 we observe interesting 



^Of course only holomorphic functions correspond to observables in the original theory. 
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Figure 4: Distribution of action density S/Q for various values of /3 at yU = on 
a 8^ lattice, using a hot and a cold start. 

crossover behaviour: both the peaked distribution bounded by i^™^^ = 6/3 and a 
decaying "tail" characteristic of small (3 distributions appear. 

To study the two possible distributions of /C™'^ further, we show in Fig. [6] 
the same results but now with /i = 0.1. In this case the hot and cold start 
yield identical distributions, since both simulations are complexified due to the 
nonzero chemical potential. The striking difference between the distributions at 
large and small /3 is still present. At large (3 the force can occasionally be large, 
making the use of an adaptive stepsize necessary. However, the typical value is 
still determined by the maximal value for real Langevin dynamics, i.e. 6/3. At 
small P this part of the distribution is completely gone and is replaced by a broad 
distribution at much larger ii'™'" values. Again at /3 = 0.5 we observe crossover 
behaviour with both features present. These results are qualitatively the same 
on larger volumes. 

Let us summarize the findings of this section. Complex Langevin dynamics 
works well at large (3 in the ordered phase. The distribution P[0\0^] in the 
complexified field space is relatively narrow in the noncompact direction and 
Langevin simulations started with hot and cold initial conditions agree. The 
drift terms do occasionally become large but the typical size is set by the maximal 
value for real Langevin evolution. At small /3, in or close to the disordered phase, 
the distribution is much wider in the 0^ direction. Typical drift terms are much 
larger, with a wide spread in the distribution. At /i = complexified dynamics 
does not reduce to real dynamics. There is a strong correlation with the phase the 
theory is in (see Fig. [2]), but not with the sign problem, since these observations 
also hold at /i = and are independent of the lattice volume. Moreover, for the 
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lattice volumes we consider the sign problem is not severe. We emphasize that 
a firm conclusion can only be drawn after all the findings presented above are 
combined consistently, while the observation of e.g. large drift terms or a large 
width by itself would clearly be insufficient. 

6 Conclusion 

We have studied the applicability of complex Langevin dynamics to simulate field 
theories with a complex action due to a finite chemical potential, in the case of 
the three-dimensional XY model. Using analytical continuation from imaginary 
chemical potential and comparison with the world line formulation we found 
that complex Langevin dynamics yields reliable results at larger P but fails when 
P < 0.5 at small chemical potential. We established that the region of failure is 
strongly correlated with the part of the phase diagram which corresponds to the 
disordered phase. We have verified that these conclusions do not depend on the 
lattice volume. Failure at small (3 values was also observed a long time ago in 
the case of SU(3) field theory in the presence of static charges [M] . 

Due to the use of an adaptive stepsize algorithm no runaways or instabilities 
have been observed. The results we found in the disordered phase are therefore 
interpreted as convergence to the wrong result. To analyze this, we have studied 
properties of the dynamics and field distributions in the complexified field space. 
For the smaller /3 values, we found that complexified dynamics does not reduce 
to real dynamics when /x = 0. Furthermore, for the system sizes and parameter 
values we used, the sign problem is not severe. We conclude therefore that the 
failure is not due to the presence of the sign problem, but rather due to an 
incorrect exploration of the complexified field space by the Langevin evolution. 
The forces appearing in the stochastic process behave very differently at large 
and small /3. Interestingly, in the crossover region at /3 ~ 0.5, the dynamics 
shows a combination of large and small /3 characteristics. It would be interesting 
to further understand this, e.g. in terms of competing (nonclassical) fixed points. 

We found that several features resemble those found in simulations of simple 
models with complex noise [l5l|60]. Our hope is therefore that a detailed study 
of simple models with complex noise can shed light on the features observed here 
with real noise. Such an investigation is currently in progress. 
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